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Abstract. Active Brownian particles (ABPs) and Run-and-Tumble par¬ 
ticles (RTPs) both self-propel at fixed speed v along a body-axis u 
that reorients either through slow angular diffusion (ABPs) or sud¬ 
den complete randomisation (RTPs). We compare the physics of these 
two model systems both at microscopic and macroscopic scales. Us¬ 
ing exact results for their steady-state distribution in the presence of 
external potentials, we show that they both admit the same effective 
equilibrium regime perturbatively that breaks down for stronger exter¬ 
nal potentials, in a model-dependent way. In the presence of collisional 
repulsions such particles slow down at high density: their propulsive 
effort is unchanged, but their average speed along u becomes v[p) < v. 
A fruitful avenue is then to construct a mean-field description in which 
particles are ghost-like and have no collisions, but swim at a variable 
speed V that is an explicit function or functional of the density p. We 
give numerical evidence that the recently shown equivalence of the fluc¬ 
tuating hydrodynamics of ABPs and RTPs in this case, which we de¬ 
tail here, extends to microscopic models of ABPs and RTPs interacting 
with repulsive forces. 


1 Introduction 

Outside the realm of exact results, and despite recent progress [1], non-equilibrium 
statistical mechanics largely remains (paraphrasing Harish-Chandra [5]) messy, elu¬ 
sive and non-rigorous. However, we may say (paraphrasing Dyson [2]) that this is 
precisely why it is such an interesting research field. While a fully general exten¬ 
sion of the well-established framework of equilibrium statistical physics remains out 
of reach, there are important classes of non-equilibrium systems for which progress 
towards a comprehensive theory seems achievable. Colloidal self-propelled particles 
(SPPs), which represent a central focus of research into active matter [3141516) . ar¬ 
guably offer such a prospect. This applies at least in some simplified limits, such as 
spherical SPPs interacting by central forces only (a restriction that excludes hydro- 
dynamic interactions). In developing theories of such systems, an important general 
question is how far macroscopic behaviour is sensitive to the detailed dynamical rules. 
For thermal systems, steady-state properties are governed by Boltzmann equilibrium, 
and hence sensitive to the energy landscape but not the dynamics used to explore it 
- so long as this dynamics obeys the principle of detailed balance (as it must do in 
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such cases). For non-equilibrium systems, such as SPPs, the question must be asked 
afresh. In this article we address this issue, by studying two classes of active particles, 

ABPs and RTPs. These have related but distinct dynamics, and have become two 
major workhorses of recent simulation and theoretical studies of active matter. We 
will find instances where their dynamical differences disappear on coarse-graining, 
to give the same macroscopic physics; and other instances where these differences 
remain important - typically when there is short-scale structure in the problem, such 
as confinement in small traps. 

In active matter, a continuous supply of energy destroys microscopic time-reversal 
symmetry (TRS) and allows phenomena to arise that are impossible in thermal equi¬ 
librium systems, where detailed balance restores TRS in the steady state. Inspired in 
large part by experiments on synthetic self-propelled colloidal particles |7I8I9I10I11I12I13I14] , 
many recent theory and simulation papers have addressed the physics of a simplified 
model comprising spherical active Brownian particles (ABPs) |15ll6ll7ll8ll9l20l21l22l28l24l25j . 
In parallel, growing experimental interest in bacterial motion |2bl27l28l29l:i(TI^Tm| 
has brought run-and-tumble particles (RTPs) to the fore as a second prototypical 
model of self-propelled particles at the colloidal scale. 

In these models, each particle feels a constant external force (vq oriented along 
its swim direction u. At low density this is equal and opposite to the local drag 
force felt by a particle moving at speed Vq- (At high density the external force may 
be partly opposed instead by conservative interparticle forces.) This is a one-body 
drag and thus the models entirely omit hydrodynamic interactions and all the other 
physics associated with the presence of an incompressible solvent. These can have 
important effects, for instance generating self-pumping states or causing ro¬ 

tation by collision |35l36j . On the contrary, in the spherical ABP model the angular 
dynamics comprises rotational diffusion with a fixed diffusivity Dr- Idealized RTPs 
differ from ABPs only in their rotational relaxation; instead of continuous diffusion, 

RTPs undergo discrete tumbling events of short duration that completely randomize 
their swimming direction; these events occur randomly in time at mean rate a and 
the mean time between two tumbles is therefore exponentially distributed Both 
models share an important simplification, which would be inadmissible for any par¬ 
ticles capable of exerting significant torques on one another (such as non-spherical 
SPPs at high density). Specifically, the angular dynamics of each particle remains 
unperturbed by interactions of any kind, and proceeds independent of the positions 
or orientations of all other particles in the system. In what follows, the terms ABP 
and RTP refer to this spherical, torque-free case unless stated otherwise. 

Another key feature of these two models is that each particle exchanges momentum 
with external driving and drag forces, not with a suspending solvent. The resulting 
non-conservation of momentum within the system greatly complicates the discussion 
of the macroscopic force balances that normally underly mechanical properties such 
as pressure I22I23I24] . Indeed, the force per unit area on a wall in general depends 
on the wall-particle interactions [SU]; this dependence does cancel out for the case of 
spherical (torque-free) ABPs moving with constant propulsive force [ID], but not if 
the propulsive force is itself density-dependent [39] (a case we return to below). 

In this paper, we present a comparative study of these simplified models of self- 
propelled particles. We first consider the steady-state distribution of non-interacting 
SPPs submitted to an external potential 14xt- Using exact results on the sedimen¬ 
tation of active particles, we show that ABPs and RTPs each admit an “effective 
temperature” regime when the Stokes speed Vg = — VI4xt/C is small compared to 
the swim speed v. We show how the effective temperature concept breaks down 


^ Other type of distributions, possibly with large tails, may also be considered. See, 
e.g., |37|38 | 
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outside this regime in model-specific ways; this issue has attracted lots of interest 
recently |l()l;UI41l4ill4;!l . We then consider the confinement of SPPs in circular or 
spherical traps, and show how the different angular dynamics between the two systems 
result in qualitatively different behaviours outside the effective equilibrium regime. 

We then turn to interacting SPPs and consider models that replace interpar¬ 
ticle collisions with a density-dependent propulsion force. Therefore we avoid any 
detailed discussion of the mechanical pressure, although it may offer an interesting 
alternative perspective on phase equilibria |24I40I44] . We explore an approximation 
scheme for torque-free SPPs with collisional interactions whereby the conservative 
forces responsible for collisional slow-down of the particles are replaced at mean-field 
level by a “programmed” slow-down, that is, an effective reduction in the propulsive 
force at high density. This path was first sketched out in |45| and pursued further in 
[17119146147148] . In a more general context it is both legitimate and interesting to con¬ 
sider, in its own right, the case with no conservative force field between particles, and 
ask about the effects of programmed slow-down among such “ghost” particles. While 
not forgetting its original motivation in the collisional case, we will often adopt this 
viewpoint here. Indeed, the concept of a programmed slow-down is well established in 
biological systems such as bacteria, where density-dependent dynamics can be effec¬ 
tively introduced through a biochemical pathway called quorum sensing [41,30149150] . 

This causes bacteria to change behaviour in response to the concentration of a short¬ 
lived chemical that they also emit. Because of the short lifetime, the response is 
then directly sensitive to the number of neighbours within a short but finite dis¬ 
tance, and hence to a local coarse-grained density. Note that this is distinct from 
chemotaxis [5n^ in which organisms swim preferentially up or down gradients of a 
long-lived chemical signal—effectively creating a long-range repulsion or attraction. 

(Chemotaxis might also arise in synthetic, self-phoretic colloids [S3].) An interesting 
study of chemotactic interactions in ABPs, showing some features resembling phase 
separation of the type discussed below, can be found in [54] . 

In bacteria such as E. coli angular relaxation is an RTP process, for which the 
consequences of density-dependent speed were first worked out in |4S|. However, some 
bacterial types, including mutant strains of E. coli known as “smooth swimmers”, 
do not tumble. In practice the genetic engineering strategies used in [30149] to give 
density-dependent slow-down primarily achieve this by altering tumble times and fre¬ 
quencies (reducing the time-averaged speed). However, were it possible to directly link 
the actual propulsion speed of bacteria to their local density in a smooth-swimming 
strain, then this would offer an unambiguous physical realization of ABPs (as opposed 
to RTFs) with a density dependent speed. As in [30149] . this could probably be made 
to happen at densities too low for collisional or crowding effects to be important. 

In summary, the study of SPPs with density-dependent propulsion speed is moti¬ 
vated (i) as an approximate representation of the collisional slow-down in SPPs with 
pairwise repulsions; (ii) as a representation of smooth-swimming or run-and-tumble 
bacteria with slow-down caused by quorum sensing or a related, non-collisional mech¬ 
anism; and (iii) as a model in its own right, with which to explore generic collective 
phenomena in active matter systems. In the last context, it offers further insights into 
the degree to which macroscopic outcomes depend on local dynamical rules, and we 
shall focus on this below in comparing the ABP and RTP cases in some detail. 

The most striking many-body phenomenon seen in this class of models is motility- 
induced phase separation, or MIPS. This is by now a well-established concept [4I45I46I47I48I55I56I57I58I5'3] 
and is the result of two effects in combination. First, the mean speed of a motile 
particle along its propulsion direction decreases with density (in contrast to the equi¬ 
librium case where the velocity statistics of a particle are separable, and hxed solely 
by equipartition). Second, particles tend to accumulate in spatial regions where they 
move slowly (a possibility ruled out, in isothermal equilibrium, by the fact that the 
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speed distribution cannot depend on position). This combination creates a positive 
feedback which is the origin of MIPS |4l45l56j . 

In what follows, after exploring in Sectionthe one-body behaviour of ABPs and 
RTPs in external potentials, we turn in Section]^ to their many-body physics, show¬ 
ing first how the large-scale dynamics with density-dependent v{p) can in each case 
be mapped on the equilibrium dynamics of passive particles with attractive forces, 
and then outlining how gradient terms affect this mapping. Sections [2| and |3| offer a 
more complete presentation of results reported in [41l46j respectively. In fact we add 
substantially to these results: our discussion of ABPs in traps and of angular distri¬ 
butions for ABPs and RTPs under gravity are both new to the current paper. In the 
many-body section we carefully compare our approach with recent work by others 
on the role of leading order nonlocality (gradient corrections) in an effective free en¬ 
ergy picture. We show further that, somewhat surprisingly, the large-scale equivalence 
between ABPs and RTPs established at the level of fluctuating hydrodynamics for 
particles with density-dependent swim speed |46j effectively extends to microscopic 
simulations of ABPs and RTPs interacting with repulsive forces. 


2 Steady-state of active particles in external potential 


Consider the classical experimental geometry, devised by Jean Perrin, in which dilute 
colloidal particles are allowed to sediment under gravity in a cuvette. Equilibrium 
statistical mechanics asserts that we can forget all details of the diffusive dynamics 
of the colloids: the steady-state probability of finding a particle at a given height z is 
a Boltzmann distribution, P{z) oc eiq){—Smgz/kT), where 6m is the mass difference 
between the particle and a corresponding volume of liquid 0 Let us now replace the 
Brownian random-walk diffusion of the colloids by an active random walk, powered by 
some microscopic non-equilibrium process. In the SPP context the latter is a persistent 
random walk of fixed speed, with either continuous (ABP) or discrete (RTP) changes 
of direction. (Its persistence length is much longer than for true Brownian motion, 
although the latter is finite in principle. More importantly, the instantaneous speed 
of true Brownian motion is not fixed, but has unbounded fluctuations given by the 
Maxwell-Boltzmann distribution.) Though not identical, clearly the equilibrium and 
non-equilibrium diffusive processes cannot be completely unrelated, at least in a limit 
where the persistence length is much smaller than the sedimentation length. Outside 
this regime, however, very different physics may be encountered. 

Previous studies of sedimenting active particles confirm this line of reasoning. 
Exact results for RTPs in one, two and three spatial dimensions have shown the 
density profile to be exponential p{z) oc exp(—Az) far away from the containing 
boundaries |41l45j . The sedimentation length reduces to = kTes/Smg, with Teg = 
v'^C/da and f the inverse mobility of the particle; this result holds in the limit where 
the sedimentation speed Vg = f~^dmg is much smaller than the swim speed v. 

This effective equilibrium regime, where the sedimentation profile is given by a 
Boltzmann weight with a temperature independent of the potential I4xt) was 
observed experimentally for ABPs in pfTilTTI and derived analytically in [^. Its 
breakdown for stronger confinement, predicted in |41j . was not observed experimen¬ 
tally so far nnEa, triggering a debate about the generic sedimen tation behaviour 
for SPPs [12]. Below, for both RTPs and ABPs (Sections 2.1|2.2 ) we derive in 2D 
the full probability distribution P(z,0) for finding a particle at height z swimming 
in direction 6. The outcome is that ABPs and RTPs both admit the same effective 


^ Note that in all this article but the sedimentation sections, we silently assume the mass 
of particles to be equal to one. 
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equilibrium regime - up to a simple parameter mapping between a and D^- However, 
beyond this regime the density profiles are dynamics-dependent, and various defini¬ 
tions of the “effective temperature”, which all coincide within the regime of small Vs, 
are found to differ beyond it. 

In Section 2.4 we turn from sedimentation (which has Vg uniform in space) to 
address the steady-state distributions of ABPs and RTFs in harmonic traps, where 
Vs{r) = —Kr. Again, ABPs and RTFs yield the same limiting effective equilibrium 
regime but lead to very different physics once the persistence length (respectively 
v/a or v/Dr) is much larger than the trap radius v/k. In this regime, the cases of 
continuous (ABP) and discrete (RTP) angular rotation give entirely different statistics 
for the density at the centre of the trap, for reasons we discuss. 


2.1 Sedimentation of 2D RTFs 


Let us consider the dynamics of a 2D RTP in the presence of an external potential 
I4xt(r) = 6mgz. The self-propulsion velocity '(;u(0) is supplemented by a sedimenta¬ 
tion velocity —Vge^ = —C~^VI4xt(r) so that the master equation reads 


P(r,0) 


-V • [(uu(6») - Vse^)P{r, 6)] - aP{r, 0) + ^ 


de'P{r, 9') 


( 1 ) 


We consider a system which is not free-falling, i.e. such that a wall at z = 0 prevents 
the particles from crossing this plane. A number of boundary conditions can be used 
to model the relevant interaction; for example, one can simply cancel the z component 
of the particles velocity m. or take into account the torque that a wall would exert 
on actual bacteria [41] . Different boundary layer structures then arise proximal to 
the wall, depending on the details of the dynamics and the boundary condition, 
which typically show strong particle accumulation in this region EEg. Beyond this 
proximal regime (typically a few run lengths in height) a stationary profile is reached 
at larger z whose form is independent of the boundary condition. 

In what follows we consider only this distal part of the profile. Its form can be 
found by assuming a factorized steady-state P{z,9) = p{z)f{9); from Eq. the 
factors must then satisfy 


P'jz) ^ a /J_ 
p(z) Vcos9 —Vs 

where A is a constant, a priori unknown. Eqs. © directly give 



p(z) = Poe 


f{e) 


1 

- ^(^cos6» - Vs)] 


( 2 ) 

(3) 


The constant A can then be fixed by the normalization condition f{9)d9 = 1, 
which yields 




IX{vs -I- u) -I- a 
A(us — u) -I- a 


= X{vs -I- u) -I- a 


(4) 


Finally, the inverse sedimentation length is given by 


2aVs 

v2 _ y2 


( 5 ) 
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Fig. 1. Sedimentation of RTFs and ABPs (a = Dr = « = !). Symbols correspond to 
simulations data and full lines to analytical predictions. The conditional probability P{6\z) 
for ABPs (top left) and RTFs (top right) is independent of 2 in the distal part of the profile. 
(For each value of Vs, several P{9\z) are constructed for different 2 in the distal region; they 
show perfect overlap.) The simulation data matches the theoretical predictions given by § 
and (lOl. Bottom left: Inverse sedimentation length A as Vs/v is varied. The equivalence 
between ABPs and RTFs breaks down as soon as one leaves the effective equilibrium regime 
(bottom right). The dashed lines correspond to the first order corrections of the two systems 
given in (12 I and (151. 


SO that 


/( 0 ) 



2tt 



2^ cos 9 + 



( 6 ) 


These exact results extend those of im where only p{z) was computed; they perfectly 
fit simulations of sedimenting RTFs, as shown in figure 


2.2 Sedimentation of 2D ABPs 

A similar path can be followed for sedimenting ABPs in 2D, starting from the master 
equation 

P(r, 6>) = - V • [{vu{e) - Vse^)P{r, 6)] + DrdeePir, 9) (7) 

Again, a separation of spatial and angular variables is confirmed nnmerically beyond 
a boundary layer proximal to the wall. The steady-state P{z, 9) is found by using the 
ansatz P{z,9) — p{z)f{9) in Eqj^ to give 

0 = -a^[(ucos6» - Vs)P{z,9)] + DrdgeP{.z,9) (8) 

which yields 

p'{z) = -\p{z) f”{9) = -^{vcos9-Vs)f{9) (9) 

As before, this yields an exponential atmosphere p{z) = po exp(—Az) where the inverse 
sedimentation length A is set by the normalisation of f{9). The equation for f{9) is 
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the Mathieu equation. Given that f{9) is even, it admits the general solution 


/W=crC(-l^,-^,2d 


( 10 ) 


where C{a^q,'d) is the cosine Mathieu function |ti4) . For fixed q, this is tr-periodic in 
-d only for a countable subset of a values, a = an{q), which satisfy a„(0) = Since 
A —>■ 0 when Vs —t 0, the solution of interest has n = 0 which implies 



( 11 ) 


Solving this numerically gives X{vs,v, Dr). (The numerics can be done using either 
a formal solver like Mathematica or by procedures based on series expansion of the 
a„’s [SS].) The resulting values of A match perfectly the results found by simulating 
sedimenting ABPs; the full distribution is then obtained by plotting the corresponding 
Mathieu function, and also agrees very well with the numerics (see Fig[^. 


2.3 Effective equilibrium regimes 


Boltzmann distributions. When Vs/v <C 1, both ABPs and RTPs admit an effec¬ 
tive equilibrium regime. This is most easily seen for RTPs where A can be readily 
expanded: 


A ~ 

Vs"^V 


^( 1 + 4 ) 


( 12 ) 


Keeping only the dominant contribution, the steady-state thus reduces to the Boltz¬ 
mann form 

p{z) = kTes = ^ (13) 

2a 

where the last equality is an effective Stokes-Einstein relation that connects the mo¬ 
bility to the diffusivity Dq = v"^/{2a) through an effective temperature Teff- 
For ABPs, one can use the expansion 


to get 




A ~ 

V3<.V 


7 4 

- Q 

128 

(14) 

7v^ 

4^2 

(15) 


Again, the system admits a Boltzmann distri 
order in Vs/v, with an effective temperature 


kTes = 


2Dr 


(16) 


Equating the mean time r = a~^ between two tumbles in RTPs and the rotational 
diffusion time r = D~^ in ABPs, the two effective equilibrium regimes match. (In 
general dimensions, the parameter mapping is from a~^ to {d— 1)D~^, which are the 
angular autocorrelation times for RTPs and ATPs respectively [M].) However, even 
the first order corrections in (12) and (15) are different in form (see Pig. [^. Match¬ 
ing the angular relaxation times gives exponential decay of the angular correlator 
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(u(i)u(O)) = exp[—t/r] in both cases, but unless Vg/v <C 1 the sedimentation length 
depends on the details of the angular dynamics. 

Note that one can also expand the angular distribution f{9) close to the effec¬ 
tive temperature regime Vg ^ v. For 2d RTFs, this gives (using the normalisation 
/<!»/(») = 1) 




cos 9 


2a2 


cos 29 + 0{k^ 


2v, 


2v 


= 1-1-cos 9 -\ -^ cos 29 + 0{[ — ) ) 


where on the last line we used that 


K = ^(1 + ^) 

Vs<.v V V 


(17) 

(18) 


(19) 


We can do the same for 2d ABPs using the Fourier expansion of characteristic 
Mathieu functions [51], to get 


-coaO + ^co»2« + 0(K=) 


2i/(9) = 1 + 

= l+?Z£cosfl+Aeos2e + 0(f^i')’) 

V 2v^ \ V / 

where on the last line we used that 


^ 2DrVs \ 

Vs-^V v‘^ 


( 20 ) 

( 21 ) 


( 22 ) 


The difference between f{9) for ABPs and RTFs can thus be seen at the level of 
the second harmonics. As for the difference between the density profiles of RTFs and 
ABPs, it appears at second order in Vg/v. Note that the first order correction to 
f{9) = l/(27r) however suffices to distinguish these distributions from equilibrium 
ones. 


Fluctuation-dissipation relations. A quantity that is frequently looked at in non¬ 
equilibrium statistical physics is the ratio between correlation and response func¬ 
tions [66] . Here we examine this for active particles in a sedimentation setup. For 
simplicity we only consider RTFs where we have explicit formulae, but conceptually 
what follows applies equally to ABPs. Let us consider a small force / applied along 
and compute the response function 


R = 


df 


f=o 


(23) 


The steady-state distribution is then 0 ( 2 ;) oc A(/) exp[—A(/) 2 :] where A(/) is obtained 
by replacing Vg ^ Vg — f in Eq. ([^. Since (z) = A(/)“^, R is given by 


R = - 


m 

xioy 


(24) 


Since the correlation function C = (z^) — (z)^ = A(0) one directly gets 


C _ 1 

h ~ ~A'(0) 


(25) 
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In the effective equilibrium regime, where A(/) = 2{vs — j this gives 


C 

h 


CV^T 


= kTeS 


(26) 


which is the usual form for a system with an effective temperature; the system satisfies 
a fluctuation dissipation theorem (as it should) governed by that temperature. 

For larger Vg , the same ratio can be computed explicitly for RTFs using the exact 
result yielding 

C — v1 

2 v'^+vl ^ 

In this regime, there is no longer an FDT at temperature Tgfj for the response to a 
force along z. Moreover, an alternative definition of the effective temperature (which 
coincides with the fluctuation-dissipation ratio in the effective equilibrium regime) 
is found by applying a force instead along e^,, and asserting a Stokes-Einstein rela¬ 
tion between the mobility and the lateral diffusivity defined via the mean-square 
displacement along x. This yields kT^ff = C,v^t just as in Eq.(26), with no de¬ 
pendence at all on Vg- For large Ug, the resulting “horizontal” effective temperature 
controls neither the density profile nor the fluctuation-dissipation in the z direction. 

In summary, outside the effective equilibrium regime, the three equally plausible 
definitions 


kTes = 


C 

R 


and fcTeff = — 


Ve^tjz) 

logP(z) 


and fcTeff = C 1™ {x^{t)/t) (28) 


yield three different effective temperatures. It follows that none is really effective, in 
the sense that if a new one is needed to describe each property measured, the concept 
of effective temperature is itself ineffective. In contrast, for both ABPs and RTFs 
within the effective equilibrium regime (arising in the limit of weak gravity) all three 
definitions of temperature coincide and the dynamics effectively reduces to an equilib¬ 
rium problem. For this reason we prefer the term “effective equilibrium” to describe 
this regime over the less explicit “effective temperature” name used earlier m- 
Interestingly it was recently shown experimentally, by measuring density profiles, 
that the addition of repulsive interactions between the particles does not immediately 
destroy the effective equilibrium regime [13]. Note that observing its breakdown ex¬ 
perimentally requires a rather large ratio of Vg/v (for Vg/v ~ 0.4, the difference in A 
is only 40% for ABFs). This may explain why this breakdown has not been reported 
experimentally so far. The more general question of effective equilibrium in active 
particles remains a subtle one. For instance, if the fixed self-propulsion speed v is re¬ 
placed by a fluctuating one, using an Ornstein-Uhlenbeck process [331331 , the effective 
equilibrium regime extends to arbitrary sedimentation speed Vg, which clearly differs 
from the RTF and ABF cases treated here. As discussed in |3T], fixed propulsion 
speed is qualitatively different from any dynamics that can sample very high speeds, 
even if these are rare (as they are in the thermal Maxwell-Boltzmann distribution). 
One reason for this is clear: if there is a fixed maximum speed, then when Vg exceeds 
this, all particles are moving only downwards and the sedimentation length must 
be strictly zero (in the absence of particle-particle interactions). Rare sampling of 
high molecular speeds is why, in thermal equilibrium, a nonzero sedimentation length 
A“^ = kT/5mg is maintained even for very large Vg = Smg/C- 


2.4 Trapping of SPPs 

Apart from RTFs in one dimension |3S] , there are no exact results available for the full 
steady-state distribution of ABFs and RTFs in a harmonic trap. While the approach 
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Fig. 2. SPPs in a harmonic trap (n = k = 1 in simulation units). As $ = kt increases, p(r) 
goes from a Gaussian centered on r ~ 0 to a distribution peaked at r ~ u/k; this holds for 
both RTPs (top left) and for ABPs (top right). As diverges, the density of RTPs in the 
bulk scales as 1/^ (bottom left), while the density of ABPs in the bulk scales exponentially 
with Dr, and the profile p(r) decreases exponentially with the distance to the trap boundary 
(bottom right). 


presented in the previous section does not apply here, because P(r, 9) is not factorized, 
some progress is still possible as we show below. The Stokes velocities of SPPs are 
now position-dependent and read Vg = —Kr/C, where r is measured from the center of 
the trap. The physics is very different from the sedimentation case since the particles 
are effectively confined to r < tt = v/ k; At r = tt the trapping and propulsive forces 
compensate exactly, hence preventing the particles from moving further outside. 

The effective equilibrium regimes correspond to r and, as expected, are the 

same for ABPs and RTPs, with Boltzmann weights for the steady-state distribution 
p{r) (see Fig. 

As the ratio <1> = kt between the persistence-length vt and the trap radius tt = vj k 
increases, the time needed for the particles to cross the trap becomes much smaller 
than the mean reorientation time. Particles thus accumulate at the outskirts of the 
trap where, on average, they point outward radially. How long they persist in this 
state depends on the dynamical details; in consequence, the density distributions of 
RTPs and ABPs differ in this regime, particularly in the central region of the trap. 
Note that, once again, this scenario differs from the fluctuating speed model of |42] 
where the steady-state in a harmonic trap is always a Gaussian, with fast enough 
particles reaching arbitrarily large distances and slow enough particles remaining in 
the center of the trap. 


RTP in a harmonic trap. RTPs pointing outward at the trap boundary undergo 
instantaneous tumble events at rate a whereafter they move across the trap along a 
certain path. Despite the trap force, and somewhat surprisingly, we show below this 
path to be a straight line segment, which fully crosses the trap for large In this 
limit, the trajectory is thus a succession of chords of the circle in 2D. The fraction of 
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time spent by a particle in the bulk then corresponds to the fraction of its run-time 
needed to reach the other side of the trap; this time is proportional to 1/^. (This 
scaling is correct although in practice the path is not followed at constant speed.) 

Beyond the scaling p{r) oc o/k in the bulk, one can furthermore compute a limiting 
shape of p{r) as a/n —t 0. Let us consider a particle at position Tq on the trap 
boundary which tumbles to make an angle 9 with the horizontal axis. Its trajectory 
until the next tumble is simply given by 

r(t) = roe”''* -f -(1 — e”''‘)(cos0, sin0) (30) 

K 


which generalizes the intuitive result r = (u/k)( 2 exp[—Kt] — 1) for rg = (u/k, 0) and 
cos 9 = —1. (The latter describes overdamped relaxation across a diameter of the 
circle to a new equilibrium position opposite the previous one, and applies in the case 
where the tumble exactly reverses the swim direction.) Eq. (30) describes a linear 
trajectory joining Tq to ^u(9), where the trap force balances again the propulsion 
force. The angular momentum L{t) = (r(t) — Tq) A rit) about the initial position rg 
satisfies L(0) = 0 and L = —kL so that Lit) = 0. The velocity r and r — rg are thus 
always parallel: the torques about rg exerted by the trap force and the propulsion 
force balance to produce a straight trajectory at an angle intermediate between u 
and rg; Only the speed |f(t)| cx exp(—Atf) is thus varying. 

Let us now consider rg = (u/k, 0). As shown in Fig. the distance to the center of 
the trap first decreases from r ~ u/k to ri„in(0) = n| co^I/k, which is reached after a 
time t = (log2)/K. The distance then increases again until (almost) reaching r v/k 
before the next tumble happens. (Since the speed |v -|- v^j vanishes as r —>■ v/k, the 
particle never reaches exactly the trap boundary.) The trajectory thus crosses twice 
any annulus of radius r > rniin(0) within the trap, where it spends a fraction of 
its duration oc a(l/|f [t+(0)]| -|- l/|r[t_(0)]|). Here, we have assumed that the total 
duration of the trajectory is l/cr and t±{9) are the times at which the particle reaches 
|r| = r. These times satisfy 


(u^ — K^r^) exp(tt;t±) = — cos9) ± y^2A2^r^-a)S~0)^r2^-rA^ 


(31) 


so that the probability of finding a particle at position r can be computed as 

p2'K—2 arccos{Kr/u) 




2 &Tccos{Kr / v) 


|i’(^i)l |r(^2)| 


d9 


(32) 


Note that p{r) formally diverges as r ^ v/k] in practice this is cut off because the 
finite duration (~ l/a) of the trajectory ensures that no particle exactly reaches the 
maximal radius r = vj k at which its net speed vanishes. The smaller a, the closer to 
the trap boundary this cut-off takes place. We have taken a out of the normalization Z 
to make the scaling of p(r) apparent; Z is then independent of a and equals Z ~ 7 r/2 . 
The bottom-left quadrant of Fig shows a very good agreement between Eq. (32) 
and simulations of RTFs for o/k < 0.1. 


ABP in a harmonic trap. Contrary to RTFs, the orientation of ABFs diffuse slowly 
so that, as 9 varies, they simply slide along the boundary of the trap, rarely visiting 
its inner region. This was used recently |67l68j to compute the distribution of ABFs 
along the boundary of small traps of arbitrary shapes, assuming that the particles 
effectively never leave the boundary of the trap. As we show below, this asumption is 
a sound one since the density of particles decreases exponentially as one moves away 
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Fig. 3. SPPs in a harmonic trap (q = Dr = 0.1, v = k = 1 in simulation units). Left 
and right plots correspond to real-space trajectories and the evolution of r{t), respectively. 
The arrows in the top left panel show the orientations of the particles. Top and bottom 
correspond to RTPs and ABPs, on similar time-scales. The dashed line in the bottom right 
panel correspond to r = n[l — Dr/{2 k)]/k.. 


from the boundary and is exponentially small in at the trap centre. Introducing 
the angle (/? between the direction of the particle and the normal to the trap e^., the 
dynamics in the (r, ip) variables become 


f = —nr + V cos <p ip = \/2I\ri{t) -sin 


ip 


(33) 


where ri{t) is a Gaussian unit white noise . Using dimensionless variables r = tDr 
f = rn/v, and r 7 (r) = riif) j, Eqs. (331 become 


f = — jj- (f — cos p) 


p = a/ 277 - —— sin p 
rDr 


(34) 


where f](T) is also a Gaussian unit white noise. The angle p thus undergoes rotational 
diffusion in an effective potential j^{l — cosp) whose amplitude diverges a.s Dr / k —t' 
0. In this limit, p oscillates around p = 0 with Gaussian fluctuations {p'^) = Dr/n. 
Consequently, r ~ 1 — Dr/(2K): the particle is almost always at the border of the trap, 
although the fluctuations of p prevent it from reaching exactly r = v/k (see Fig. 
The density of particles in the bulk of the trap is thus much smaller than in the RTF 
case since in a time t ~ 1/Dr, the particles are exponentially unlikely to escape the 
trap boundary while RTPs would typically cross the trap in times f ~ 1/a. (Recall 
that these two times are interchangeable in the effective equilibrium regime.) The 
occurrence of an effective potential, preventing particles from diffusing in angle away 
from the outward normal direction, at first seems to violate the fact that the dynamics 
of ABPs is torque-free. Note however that there is no real potential for the angle d; 
the effective potential for p arises not because there is a torque, but because the 
translational motion of a particle along the trap boundary naturally leads it towards 
a location where the outward normal is parallel to the current orientation, hence 
leading to v? ~ 0. This effect is also present for RTPs but the discrete angular dynamics 
effectively allows the RTP to immediately escape the reorientation potential. 

To reach the bulk part of the trap, Eq. (34) tells us that p needs to climb the effec¬ 


tive potential to turn and face its destination. In principle, there are many stochastic 
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paths leading from the boundary to some inner target point at distance r from the 
center. As k/D j. —>■ 0, however, the transition probability is dominated by the most 
probable path leading to such a point [M] (an instanton). The probability of this path 
will be given by exp nAE/Dr, where AE(r /tt) is a geometric function depending on 
the exact location of the target. Since reaching the center of the trap requires a fluc¬ 
tuation of (p much larger than any required to reach another point close to the trap 
boundary, the energy barrier increases as the distance of a target from the boundary 
increases. In practice, a linear function AE(r Itt) fits very well the effective energy 
barrier: 

P^'^) - [C‘^[r/rT - 1)] (35) 

where C is a dimensionless constant. The density in the center of the trap indeed 
decreases much faster than for RTFs, as <?exp(—C^). Somewhat surprisingly, this 
simple argument gives reasonably good agreement, with a geometric factor C whose 
fit yields C ~ 2.15, with the bulk density measured in numerical simulations (See 

Fig.§. 


2.5 Effective equilibrium for generic potentials 

In this section we show that for generic potentials I4xt) sufficient conditions to observe 
an effective equilibrium regime are given by 


1^ ^VI4xt(r)| < v\ and Ml4xt(r)| < a, Dr{d- 1) 


(36) 


To show this, let us consider a system for which (36) holds. The fluctuating hydro¬ 
dynamics for ABPs and RTFs in the presence of external drifts derived in [3^ reads: 


p = -f V(np) -f V • [r'VFextp] + V(AVp) 
Apa = -Va{vp/Q) + V6(e-lVfaI4xtFa) + 0{V^) 


(37) 

(38) 


where fl is the surface of the unit sphere in d dimensions and A equals a for RTFs and 
Dj.{d—1) for ABFs. This description is valid at times t ^ A~^ and for fields p, p whose 
gradients are small on the scale of the persistence length of SPPs, |Vp|/p, |Vp|/|p| ^ 
A/v; its derivation in the isotropic case is the subject of section]^ whose notation 
differs slightly through introduction of a quantity p = p/ 

Expanding the second term of the r.h.s of Eq. (38) then yields 


(A - C"AVe^t)Pa = -Vaivp/n) + T'V^Fext ' ^bPa + 0(V^ 


(39) 


be neglected and (39) simplifies into 


Equation ( [39| ) thus shows that pa — 0{V). Using ( pI6| , the Laplacian of Fext can then 

(40) 


Pa = -jVa(p/n)+0(V^) 


The dynamics of p is then given by 


p = V 


(A + ^)Vp + r^VFextP 


(41) 


which is nothing but the equilibrium Fokker-Planck equation of a colloid of mobility 
^ subject to an external potential Fixt, at a temperature 


fcTeff = C[A + 


(42) 
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The steady-state solution of this equation is the Boltzmann weight p{x) oc exp[—/3effKxt], 
which satisfies 

V Vp 

A p 


V r^VKxt 
A Df +4 


< 


r'vKxt 


(43) 


which is indeed small, thanks to (36). 

The two conditions (36) are enough to show the long time dynamics of the SPP 
to amount to the equilibrium dynamics of a colloid subject to an external potential 
14xf Physically, they require that the Stokes speed Vg = —WV be much smaller that 
the self-propulsion speed v of the SPP, but also that its variations over a run-length 
Vvs-v/A be much smaller than v. 


3 Fluctuating Hydrodynamics of ABPs and RTPs 


Following [46], we address a model using a general angular relaxation dynamics, com¬ 
bining smooth angular diffusion with diffusivity Dj. with pointwise random reorienta¬ 
tions (tumbles) at rate a. The model has ABPs as one limit and RTPs as another. An 
intrinsic translational diffusivity Dt is also allowed for; the swim speed v is allowed to 
be position-dependent in the first part of the derivation, w(r), and we will then turn 
to density-dependent v{p). For a single particle, the probability density tj}{r,u,t) of 
finding the particle at position r moving in direction u obeys exactly (in d = 2,3) 

'0(r, u) = -V • [uuV'(r, u)] -b Vu • [Dr^uil^ir, u)] 

a f (44) 

-b V • (AV'(/'(r, u)) - a'!/'(r, u)-b — / ^(r,u')dl7' 

where Vu is the rotational gradient acting on u and the integral is over the unit sphere 
|u'| = 1 of area 17. The first term on the right is the divergence of the advective current 
resulting from self propulsion and the last two terms are loss and gain due to tumbling 
out of and into the direction u. The second and third terms account for rotational 
and translational diffusion. 


3.1 Coarse-graining procedure 


The first part of the derivation below consists in using a moment expansion to show 
that at large time and space scales, the dynamics (44) amount to a Langevin equa¬ 
tion, whose drift and diffusion terms can be fully characterised. Using ltd calculus, 
we will then start from this single-particle Langevin dynamics to derive a collective 
Langevin dynamics of the density field of N non-interacting active particles. Since all 
the computations will be done allowing for dependence of the microscopic parameters 
V, a, Dr on the spatial position r of the active particle, our computation applies to the 
case where these dependencies occur through the density field p{r). This derivation 
will thus provide the stochastic dynamics for the density field of N interacting active 
particles. Were we to consider purely ABPs, one could directly use Ito calculus to 
construct the dynamics of the density field, as was done in 155170171 j . In the general 
case, it is the presence of the tumbles that makes it necessary to take a first diffusive 
limit. Note that in lattice models, field theoretic methods can be used directly that 
bypass this limitation and give the same final result m- 
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3.1.1 Moment expansion 
We first decompose ip as 

^p{r, u,t) = (f + p.u + Q ; (uu - I/d) + 


(45) 


Here, (/j, p, Q are functions of which parameterize the zeroth, first and second 
angular (d = 2) or spherical (d = 3) harmonic components of '0, while O projects onto 
the higher harmonics. Note that in 2 dimensions, the angular harmonics are Fourier 
sine and cosine series in 9; clearly 

• a , n J ^ \ Qll ~ ^22 , *5l2 + Q 2 I . 

p.u = Pi smd+p 2 cosd and ii:(UiUj —-2-cos2dH- - -sm2d 

span the subspaces generated by the first and second harmonics, respectively. The 
same holds in d = 3, where the components of any unit vector Ui are linear combina¬ 
tions of y/", m = —1, 0,1 whereas the components of the traceless tensor UiUj — 6ij/3 
are linear combinations of Y^, m = —2... + 2. 

We now proceed order by order in the harmonics. Integrating Eq. (44) over u gives 

(46) 


1 . 


p = -^V(up) -b V(AVp) 


whereas multiplying eq. (44) by u and then integrating over u gives 
Pa = -Va(np) - {Dr{d- 1) -|- a)pa -f V(AVpa) “ 


Finally, multiplying eq. (44) by uu — I/d and integrating over u yields 
d-b2 


Qab — 


^-BabcdScVPd - {2dDr + a)Qab + V • [DtVQab] - Y'cVXabc ( 48 ) 


where Babcd = {dacha + dadhc - 2Sab6cd/d)/{d + 2) and Xabc, which comes from 
highe r order harmonics, will not play any role in the following. The derivations of 
Eqs. (46p8) are detailed in Appendix [A} 


3.1.2 Diffusion-drift equations 


So far, beyond the assumed isotropy of u(r), Z?r,t(r) and Q!(r), no approximation has 
been made; Eqs. ( 46p8 ) are exact results for the time evolution of the zeroth, first 
and second harmonicsm ip{r, u, t). We then note that Eq. ( [4^ is a mass conservation 
equation: 

Lf) = —V • J with J = -vp — DtVip (49) 


with a current J which involves the first harmonic p. We will now use a large time 
and space scale limit to obtain J as a function of p. 

We first note that p is locally conserved and is thus a slow mode: the relaxation 
of a density perturbation on a scale £ occurs in a time that diverges as (. —)■ oo. 
On the contrary, p and Q are fast modes, whose relaxation rates are given by tp = 
a -b Dr{d — 1) and Iq = a + 2dDr. For times greater than fp.q, one thus assumes 


p = Q = 6>[i/] = 0. Eqs. (47) and (48) then give the ordinary differential equations 
that p and Q satisfy quasi-statically as p evolves on time-scales much larger than 
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tp.Q. By itself, this creates a non-local constitutive relation between J and tp which 
still involves all harmonics. 

Next we explicitly carry out a gradient expansion, yielding for Q: 


Qab — 


1 


2dDr 


d + 2 


BabccNcivpd) + Vc(Xabc) + C>(V 


from which the quasi-stationary p then follows via (|47| as 

p = -. 


1 


{d — l)Dr -b a 


V(t;(^) + 0(V") 


(50) 


(51) 


At this order in the gradient expansion, closure is achieved without needing further 
information on harmonics beyond the first; the current J is given by 


J = - 


did — l)Dr + da 


V(vip) - DtVip 


(52) 


At this order, Eq.(49l corresponds to a diffusion-drift approximation of the micro- 

(53) 


scopic master equation (44), given by: 

ip = —V • [—D'S/p + V(p] 


where the diffusivity and drift velocity obey 

—uVu 

did — l)Dr da~^ * ^ did—l)Dr + da 


(54) 


At this level of description, though not for the exact results that preceded it, RTFs and 
ABPs are seen to be equivalent. That is, in (53), the tumble rate a and rotational 
diffusivity Dt enter only through the combination [d — 1)D^ + a, so that the two 
types of angular relaxation are fully interchangeable. Conversely, given a sequence of 
snapshots showing the large-scale evolution of the local densities of SPPs, one cannot 
determine whether the system is composed of ABPs or RTFs. Note that “large scale” 
here excludes some of the trap problems considered in Sectionj^where the confinement 
length is smaller than the persistence length of the self-propelled motion. 


3.1.3 Many-body physics 


Eqs. (49) and (53) give the evolution for the probability density of one particle at 
Following [45j we can now consider an assembly of particles whose 
and Dt depend on position both directly and indirectly. 


diffusion-drift leve 
motility parameters v, a, D. 


through a functional dependence on the microscopic density p(r, t) ~ J2f_i 
Since we are considering a multiplicative noise whose variance is a functional of the 
density field p, the construction of the Langevin equation associated to p is technically 
more involved than in the additive case m . First, we go from the Fokker-Planck equa¬ 
tion (531 to the equivalent Itb-Langevin equation for an individual particle position 

<5 


,(t) = V + Vii?(r^,[p]) + (Vr 




)Divti,[p\) + ^^'n 


(55) 


where Vil?(r^,[p]) represents a gradient with respect to the first argument of D 
and not with respect to its implicit dependence on through the field p(r') = 
5(r' — r„). On the contrary, the operator Vr^ ^ applies to the functional 
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dependence of D on the field p{r') and the subscript shows Vr^ to apply to 
dependencies on and not on r' or that may appear in D. (For more technical 
details, see Appendixh^) The two terms Vil?(r^, [p]) + V [p]) form the 
so-called “spurious d^” term due to the Ito time discretisation. (This term is, of 
course, not spurious, but necessary once one uses the Ito convention to define that 
discretization [7371].) The atypical form of this spurious drift is due to D depending 
on both explicitly, which explains the first term, and implicitly—through p(r')— 
which explains the second one. In Appen dix [B] we show that the density p(r, t) then 
obeys the many-body Langevin equation [45] 


p = -V. (^V[p]p - D[p\Vp + (Vr^)i?(r, [p]) + {2Dpf/^^ 


(56) 


with white noise {Ai{r, t)Aj(r', t')) = 8ij5{v — r')5{t — t'). The gradients in ( [^ have 
no subscript ‘I’ since there are no more ambiguities at this stage: the derivatives are 
taken with respect to the spatial coordinate r, which does not enter p(r'), and not 
with respect to the position of one of the N particles as in (55). Note that when the 


diffusivity of particle p does not include ‘self-interaction’, i.e. D — I?(r^, [p—5(r —r^)] 
or is the result of a convolution between p(r') and a symmetric kernel K{v), i.e. D = 
f dr'A:(r^ - r')p(r'), then [p]) = -VA:(0) vanishes (see Appendixj^ 

for details) and one recovers the more standard equation 


p = -V. (V[p]p - D[p]Vp + (2Z?p)1/2a) 


(57) 


The functionals u[p], a[p] and 7It,r-[p] in ( [54] ) then define for the interacting particle 
system the many-body drift velocity and diffusivity V[p] and D[p\ for use in Eq. (56). 

Starting from the many-body Itb-Langevin equation ( [5^ , one can then derive 
a functional Fokker-Planck equation for the evolution of the probability density 
P[p{r),t] of finding the system with a density field p(r) at time t: 


P[p] = 



pV — HVp — Dp 



P[P] 


(58) 


The technical details, which show the importance of the atypical spurious drift, are 
detailed in appendix [C| 


3.1.4 Connection to large-deviation functionals; role of noise 


A crucial observation, first made for RTFs in [33 ; is that ( [5^ reduces, under specific 
conditions, to a description of passive Brownian particlesfPBPs) with a specified 
free energy functional fiP[p] = /3J^ex[p] + / p(lnp — l)dx. (In what follows we use 
thermal units in which ,0 = 1.) This is best seen by noting that the Fokker-Planck 
equation (58) admits flux-free solution 


pV — HVp — Dp 



P[p\ = 0 


whenever there exists a functional Pe-x[p] which satisfies the condition 


(59) 


V([p],r)/£>([p],r) = -Vr((5JexM/i5p(i-)) 


(60) 
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Indeed, in such a case P[p] = exp[—/3J^[yo]] satisfies (f^. As was noted in |71l47j . 
s[p] = ISPjV, with V the volume of the system, is the large deviation function of the 
density profile p 


s[p] 


lim ^ 

V —>-oo V 


logP[p] 


(61) 


It is at first sight very surprising that one can compute this object for any kind 
of interacting out-of-equilibrium system. However its computability depends on (60) 
holding true, and indeed this is not true in any general manner. It is however satisfied 
to leading (zeroth) order in a gradient expansion of Tex- If (60) holds, the system is 
completely equivalent to an equilibrium system of passive particles with excess chem¬ 
ical potential gradient (57ex/^p(r)- This equivalence holds not only at the level of the 


steady-state but also at the level of the dynamics (58). This means that, at this macro¬ 
scopic level, the dynamics is effectively an equilibrium one and, for instance, satisfy 
the Onsager-Machlup symmetry between excursion and relaxation |69j . However this 
does not mean that the phenomenology of the system reduces to an equilibrium one. 
For instance, the slowing down of active particles at high density causes an effective 
many-body attraction between the equivalent passive particles, characterized by a 
J^ex whose local part has negative curvature in p. This kinetic slowdown can trigger 
a phase separation mechanism, which we detail in the following section, which would 
be impossible at equilibrium. (Making the viscous drag on each particle an increasing 
function of local density would generate a similar kinetic slowdown in an equilibrium 
system, but the Boltzmann distribution would be insensitive to this effect.) 

If the condition (60) is not met, however, the dynamics is not equivalent to PBPs 
with conservative interactions and is therefore “irreducibly” active. As we will show 
later, the equivalence established above will only hold for homogeneous systems and, 
when the kinetic slow-down triggers a phase-separation, higher-order gradient terms 
comes into play at the large-deviation level, and break the mapping to equilibrium. 


3.1.5 Mapping to equilibrium 


The simplest case is where Dt = 0. Here the left hand side of (60) is — Vlnu[yo] and 
we then require 5Pex[p]/5p{r) = ln(u([p];r). The simplest first approach is to assume 
that the functional dependence of swim speed on density is strictly local, so that 
u([p];r) = v{p{v)). We then have = f fex(p(r))dr where f^x = fglnv(A)dA. This 
structure in the free energy is equivalent to having a passive system whose chemical 
potential obeys p = ln(/9v); the mean particle current then obeys 


J = —pDW In(pti) 


(62) 


Noting that in our chosen units (/3 = 1) the mobility coincides with the diffusivity D, 
this is the expected form for the stated chemical potential. When Dta is a nonzero 
constant, the result for generalizes to 

fex = 7; [ ln[u(A)^ -f dDta]dX (63) 

^ Jo 

This case is explored in [56] . 

To approximate the excess free energy by a local function is a widely used ap¬ 
proximation in equilibrium systems which makes sense at (Landau) mean-field level 
where fluctuations are neglected. (If the excess free energy were genuinely local, short- 
length scale fluctuations would be out of control, so to go beyond mean field theory 
requires treatment of nonlocality even if that is weak). The chosen form (with Dt = 0) 
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leads to a spinodal instability whenever dv/dp < —vjp [45156] . The system is then 
equivalent to an equilibrium system undergoing liquid-gas phase separation due to 
attractive forces and the spinodal instability is complemented by binodals, which can 
be characterised analytically in some cases (See [56] for a recent review). 

For ABPs with collisions, one finds empirically that ^(p) ~ uo(l — p/p*) with p* 
a near-close-packed density beyond which self-propulsion is effectively arrested. The 
system is then spinodally unstable for p > p*/2 [15117119] . Stability is restored by 
packing constraints at high enough density, but in the present approach (collisions 
replaced by an effective density-dependent propulsion speed), empirical corrections 
to Tex are needed to account for this HZ]. An alternative is to impose a sharp cutoff 
so that / = -boo for p > p*; this is used in m- 

The presence of a finite Dt, either due to Brownian motion or resulting from the 
random collision of particles (in which case Dt will depend on p), alters the shape 
of the spinodal curves which instead obey 1 + pv'/v < —dDt/{v^T). This causes 
the spinodals to meet at a critical point; while in qualitative agreement with the 
shape of the phase diagram reported for ABPs with repulsions (see, e.g. mnm), 
this is somewhat accidental since the actual role of Dt in the simulated dynamics is 
negligible. The observed critical point is better viewed in terms of the competition 
between slowing down at high density (promoting motility-induced phase separation) 
and the buildup of particle mechanical pressure [40] . The latter stems primarily from 
the pair interaction, rather than the small ideal gas part which is proportional to Dt 
[24l¥i] . Numerically, the spinodals are hard to locate precisely, whereas binodals can 
be located by looking directly at the coexisting densities in phase-separated states. 
The point at which phase separate is seen kinetically in ABP simulations often lies 
in between the spinodal and binodal; it depends in general on nucleation rates and is 
subject to large finite-size corrections. Very large scale simulations are thus required 
for accurate phase diagram determination mm- 


3.1.6 Beyond the local approximation 

When a system phase separates, large gradients develop and the gradient expansion 
cannot be truncated at lowest order: the local approximation to T no longer yields a 
good approximation to the large deviation functional. In such a case, one has to look 
for higher order gradient terms, which stem from two different mechanisms. First, 
for real particles, the interactions between particles are never perfectly local: even 
for hard-core repulsions, the particle size defines a finite interaction range. Second, 
we have dropped higher angular harmonic contributions by appeal to a gradient 
approximation. Retaining higher order gradients thus requires two different type of 
terms whose impact on motility-induced phase separation we now discuss. 


Non-local v{p). The local form of the free energy can be used to predict the binodal 
densities for phase coexistence via the the common tangent construction on / = 
fex + p(lnp — 1) [ISj. In this construction one equates the chemical potential df/dp 
and the ‘thermodynamic’ pressure pdf /dp — f in the two phases; this means that a 
single line can be drawn on a plot of f{p) that is tangent at the binodal densities and 
lies below / everywhere else. However, there is a hidden pitfall here: this construction 
implicitly assumes that, whatever the nonlocal terms are, these continue to obey (60). 
Whenever they do, a free energy still exists and, so long as it does, the binodal 
densities do not depend on the precise form of the nonlocal terms. Conversely, if 
a free energy does not exist, then the binodal conditions can depend explicitly on 
the nonlocal terms. This is explored in some detail, within the simplification of a (j)'^ 
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effective field theory, in [48]. For ABPs, numerical evidence of the breakdown of the 
common tangent construction was reported in mm- 

The lowest order nonlocal theory can be constructed by assuming that the func¬ 
tional dependence of swim speed on density takes the form [ 3 ^ 

v[p] = v{p) ~ v{p -I- ~ v{p) -I- u'(p) 7 ^V^p (64) 


This represents a quasi-local dependence on a quantity p that samples the local den¬ 
sity isotropically in a region of size 7 . For programmed slowing-down (e.g., quorum¬ 
sensing) the simplest assumption is that 7 is a density-independent constant. However 
for collisional slowing down in ABPs, a better estimate is 7 ~ 7 ou(p)/[(d— l)Dr -ba] 
which, with 70 of order unity, is the distance travelled during one angular relaxation 
time. Substituting the nonlocal form for v in (621 gives to leading order 


J = —pDVp 

(65) 

p. = ln(p'i;) — — 

( 66 ) 


Only if is a constant, independent of density p, does this form of chemical 

potential support the existence of a free energy. If this combination is constant - such 
as for the case where 7 is constant and v{p) ~ exp[—ap] [45] - then the common 
tangent can legitimately be used to predict the binodal densities. In all other cases 
it cannot be relied upon fT7l4^ . Note that in practice, only quantitative differences 
with the local theory have been noticed when simulating models with such non-local 
u{[p]) |17l48j : the phase separation still occurs, the coarsening law is not much affected, 
and one only observes quantitative shifts of the binodals. Importantly, the concept of 
a binodal is maintained: the densities of coexisting phases do not change as the global 
average density in the system is varied between the two coexistence values. 


Other gradient terms. The second source of higher order gradient terms is the gradi¬ 
ent expansion used to close the spherical harmonics expansion in Section [3.1.2| Were 
we to pursue this gradient expansion to higher orders, we would obtain a set of o rdi- 
nary differential equations for p, Q, etc. rather than the simple algebraic relation (511 
between p, cp and Vp. Solving these equations in terms of tp and reinjecting into J 
would then not lead to a simple Fokker-Planck equation, from which we would not 
be able to derive a microscopic Langevin equation like (55), which is the starting 
point of our approach. The equivalence between ABPs and RTPs at this higher order 
in gradients is thus questionable. This echoes the fact that the difference between 
the two models becomes more important at short length scales. As we will show in 
Section [3.1.71 however, these differences hardly impact the phase diagram. 

For ABPs, since one directly starts at the microscopic level with coupled Langevin 
equations, one can bypass the construction of the microscopic Langevin equation (551 
and directly use Ito calculus to obtain a stochastic equation for the probability density 
p(r,u) of finding particles at r going in the direction u [55170171] . One then has to 
project this equation onto successive harmonics and use an appropriate truncation. 
Again, the equation (56) amounts to second order gradient expansion. It is then easier 
to pursue this development to higher orders, which would yield coupled stochastic 
equations for p, p, Q, etc.. As far as we are aware, this has not be carried out in the 
literature for ABPs or RTPs, but progress along this path can be found for aligning 
nematic particles m- 
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For ABPs, a similar path has been followed at mean-field level, albeit without a 
proper derivation of the noise terms, to arrive at |lhl2hlK!^| 


p = DAp — V • [up] 
p = -Drip - ^V[up] -1- DAp 


(67) 


One can then always complete these equation by adding ad-hoc Gaussian noise m 
or by assuming that, were the coupled equations ( |67| ) to be written in a gradient form 
(implying the existence of a free energy functional) the noise terms would be whatever 
is required to give an equilibrium like-dynamics. As we have shown in Section [3.1.4[ 
however, the derivation of the noise terms is a crucial step to determine whether 
or not the system really maps onto an equilibrium one. Indeed, the integrability 
criterion (601 is derived by requiring a flux-free steady-state in the Fokker-Planck 


equation. This condition crucially relies on the fact that the variance of the noise 
obeys a fluctuation-dissipation relation with the mobility. 

Nevertheless, since the derivations of these noise terms are particularly difficult, a 
first strategy can be to look at whether the equations (67) have the gradient structure 
of an equilibrium model. While analyzing the linear stability of (67) for a linearly 


decreasing v{p) reveals a spinodal decomposition scenario very similar to the one 
described in section (3.1.5), it was explicitly shown in [59] that these equations are 
not of gradient form. However, expanding close to the linear instability, the authors 
of [59] were able to map these equations onto a Cahn-Hilliard model, but with an 
effective free energy which depends on where in the phase diagram the expansion is 
carried out. This revealed a rather surprising feature: the addition of the zip terms 
breaks the global mapping to equilibrium but seems to preserve a local one. This 
however fundamentally changes the structure of the phase diagram: there is nothing 
to guarantee that the binodals stay in fixed positions as one moves along the ‘tie line’ 
between them (and without this property the concept of a binodal is inapplicable). 
Moreover, according to this calculation, a large part of the transition line is changed 
from a first order to a second order transition, in contrast to the isolated critical point 
usually seen in liquid-gas type phase-separation. 


To test the predictions of this approach, we simulated directly the equations (67) 
using spectral-methods and semi-implicit time-stepping. We chose a form of v{p) 
which does not lead to v{p) < 0 at any density and thus does not require an ad-hoc 
cut-off at the level of the effective free energies: 


u(p)2 =vl + {vl - ?;g)(l - e 


( 68 ) 


For this choice, we can compute analytically spinodals and binodals predicted by the 
local theory [^ and hence evaluate precisely the impact of the Z\p term. The sce¬ 
nario revealed by our simulations does not show any of the surprising phenomenology 
predicted by the methods of jSSj: the binodals are quantitatively shifted by this new 
gradient term, in qualitative agreement with the effect of gradient terms stemming 
from a non-local v{p) [48], but the binodals do not vary along the tie-line in the co¬ 
existence region (where the lever rule still applies). Also, the transition line seems to 
remain first-order on approach to the critical point, exhibiting the familiar hysteresis 
curves (see fig. |^. 

In principle, this outcome could depend on our choice of v{p); note that direct 
comparison with [59] is not possible since the theory developed there relies on the 
addition of confining terms at the level of the free energy that do not exist at the level 
of the PDF, Eq. (67). However this interpretation seems unlikely given the growing 


breadth of literature supporting the existence of well-dehned binodals in MIPS. An 
alternative possibility is that the new and unusual phenomenology predicted from the 
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Fig. 4. The phase diagram (left) predicted by the theory without the zip term (red and 
green lines) are only qualitatively affected by the Zip term (blue). Hysteresis loops (right), 
showing the first order nature of the transition, can be seen on approach to the critical point 
by ramping up (solid lines) or down (connected symbols). D = 0.25, vi = 0.25, (f) = A, t = 
l/Dr = 1 . 



Fig. 5. Coarsening dynamics of RTFs (top) and ABPs (bottom) interacting via the density- 
dependent swim-speed (681 with vo = 15, vi = 0.25, ip = 4, t = 1, Dt = 0.25. Simulated for 
N = 120 000 particles in a box of side L = 100. 


approach of |59| is an artefact of a quasi-linear, noiseless treatement of what is in fact 
a noisy nonlinear transition. Further study is needed to resolve this issue, for inst anc e 
by carrying out the analytical procedure in for the particular v(p) chosen in (68). 


3.1.7 Comparison between microscopic simulations of RTFs and ABPs 

Since gradient terms can alter the equivalence between ABPs and RTFs, we carried 
out extensive microscopic simulations of MIPS for both models in two dimensions for 
JV particles in an L x L square domain. We first compared the two models in the 
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Fig. 6. Left: The local density distributions P{p) of ABPs and RTFs interacting via the 
density-dependent swim speed Eq. (681 with vq = 15, vi = 0.25, ip = A, t = 1, Dt = 0.25. 
The two systems overlap very well, showing that not only the mean coexistence densities 
in the two models but also their fluctuations are very similar. As expected, changing the 
mean density po does not change the coexisting densities. Right: Phase diagrams of RTFs 
and ABPs interacting via WCA potential. The boundary lines are obtained by computing 
the binodals as a function of v in the phase coexistence region for various global densities 
po = N/L^, to check that the binodals are indeed independent of po. Again, a surprisingly 
good overlap between the two systems is observed. 


case of a density-dependent swim speed v{p) given by ( 681 . The coarsening of both 
phase separating systems is shown in figurej^and shows a strikingly similar dynamics. 
Beyond the apparent similarity of the two dynamics, one can compute the steady- 
state density distributions P{p) of the two systems in the phase separated region, 
which shows that both the mean densities and their fluctuations in each phase are 
also very similar (see Fig. |^. Note that as expected, the coexisting densities along 
the tie-line in the coexistence region do not depend on the mean density po = N/Li^. 

We then simulated both ABPs and RTFs in the case where slowing down is caused 
by repulsive interactions (represented as a WCA potential as in |15I16I17] '). We know 
that the v{p) theory correctly captures the slowdown in these models but that it 
lacks the direct interparticle forces which are responsible for saturating density in 
the high-density phase [JD] . Therefore, the dynamical equivalence of ABPs and RTFs 
at large scales might break down for the repulsive case where much of the physics 
depends on short-distance collisional events taking place below the coarse-graining 
length scale of the fluctuating hydrodynamics theory. Suprisingly however, the phase 
diagrams of repulsive ABPs and RTFs collapse onto each other upon the usual rescal¬ 
ing a = (d— l)Dr (see Fig. |^. Thus, although we have shown that the details of the 
angular dynamics are important for SPPs confined to traps of size comparable to the 
run-length or smaller, the presence of short-length scale physics in the mechanism 
for MIPS does not create a similar dependence on details, as judged either by the 
qualitative kinetic observations or a quantitative study of phase diagrams. 


4 Conclusions 

Motivated by a wish to understand for active systems the relation between micro¬ 
scopic dynamics and macroscopic behaviour, we have in this paper compared in detail 
two distinct but related classes of self-propelled particles: Active Brownian particles 
and run-and-tumble particles. These differ only in the details of their angular relax¬ 
ation (continuous versus discrete). We considered first one-body-problems involving 
non-interacting particles in either uniform fields (e.g., gravity) or isotropic harmonic 
traps. In both cases, for weak enough forces there is an “effective equilibrium” regime 
in which the characteristic length scale set by the balance of external and propulsive 
forces is large compared to the run length, defined as the distance a free particle 
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can travel during its orientational relaxation time. In this regime, the two types of 
dynamics (ABP and RTP) are equivalent modulo a simple mapping between param¬ 
eters, which corresponds to equating the angular relaxation time in the two cases. 
However, for stronger confinement the corrections to this picture are specific to the 
chosen dynamics. In the case of harmonic traps, this difference is particularly striking: 
in the regime where the confinement length is large compared to the run length, there 
is an exponentially strong suppression of the density of ABPs at the trap centre when 
compared to RTPs of the same rotational relaxation time. This can be attributed to 
a subtle but physically simple mechanism whereby, in the strong trapping regime, 
particles tend to remain a long time at the outer edge of the trap, where external 
and propulsive forces balance. For RTPs this balance is suddenly destroyed with each 
tumble event, whereas for ABPs, which update their orientations gradually, there is a 
tendency to drift around the perimeter so that the balance is continually maintained, 
causing only an exponentially small escape rate from the surface region. We also con¬ 
sidered the effect of the microscopic rotational dynamics on MIPs, or motility-induced 
phase separation. For this process we gave a fuller presentation of results summarized 
in our earlier papers concerning the derivation of collective equations for the particle 
density field, from which criteria for MIPS are easily found within a limiting approx¬ 
imation corresponding to neglect of gradient terms in the large deviation functional. 
This criterion is again insensitive to the choice of rotational dynamics. It is explicitly 
derived from a model where particles interact through a programmed dependence of 
their propulsion speed on the local density, but also describes, to reasonable accu¬ 
racy, the case where slowing down is caused instead by collisions. We then considered 
the role of gradient corrections, arising either from nonlocality in the dependence of 
speed on density, or as additional terms in a combined expansion in spatial gradients 
and angular harmonics of the local distribution of particle orientations. Although 
the latter could in principle cause qualitative shifts in phase behaviour, in common 
with previous studies we do not find compelling evidence for anything more than a 
quantitative shift numerically. Specifically, although the “equilibrium” conditions for 
phase coexistence (those derivable from the local part of the large deviation func¬ 
tional) are violated, the system still exhibits conventional binodals whose defining 
property is that the density of coexisting phases is independent of the intermediate 
global density of the system. For the case where MIPS is caused by collisional rather 
than programmed speed reduction, we find, somewhat surprisingly, that there is once 
again almost no difference in phase behaviour between ABPs and RTPs with matched 
angular relaxation time. This is despite the fact that both the collisional dynamics, 
and the presence of sharp interfaces in the system, involve length scales below those 
at which the equivalence of the two models can be formally established by systematic 
coarse graining. 
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A Moment-expansion 


We start from the master equation on the probability density tp{r, u, t): 




— V • [uui/:] -I- + V • [DtVf/'] — Oilp + 


a 

7 ? 


/ 


df2il’(u') 


(69) 
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with Z\u the rotational part of the Laplacian acting on u. We then use the decompo¬ 
sition (45) of ip, which we recall here for clarity 


V’(r, u, t) = (/3-|-p-u-|-Q: (uu -1/d) + O/i/] 


(70) 


In the following we will use the notation Mab = UaUb — Sab/d and the convention 
that repeated indices are implictly summed upon: figi = fi9i- Because Q enters 
ip solely through the combination 'Y/i it can always be chosen traceless 

symmetric. (Indeed, replacing Q hy (Q + Q^)/2 ov Q — lTr{Q)/d does not change 
this sum.) 

We will use the standard notation for the scalar product on the sphere 

{f.9)=j^^f{'u)g{u) (71) 

for which spherical (and angular) harmonics are orthogonal. To obtain equations for 
p, Q, we take the scalar product of equation ( [^ with 1, u and M. To do so, we 
compute 


(1) V^) — (u, '0) — i^abi '0) — ^ BabcdQcd — ^ ^ 2 


which rely on 


/ ■ 


{uiUj,UkUe) = 


n 


d{d -I- 2) 




Bijki — ^ i^ij: ^kl) — d 2 ^ii^jk ^^ij^kl) 


(72) 


(73) 


Projecting (69) on 1 then yields for 0 


1 . 


0 = -^V[up] -f V • [AV(/3] 


while projecting (69) on u yields for p 

^ = -(u,v[uu0]) -bT>r(u, Z\u0) + V • [AV(u, 0)] - Q:(u, 0) 


(74) 


(75) 


The last two terms are easy to compute since (u, 0) = J7p/(i. Then, spherical and 
angular harmonics of order £ are eigenvectors of Z\u with eigenvalues —£{i + 1) and 
—£^, respectively, and the projection on u selects the p.u term so that 


T)r(u, Z\„0) = -Dr{d- l)f2^ 

d 


(76) 


The first term of the r.h.s. of (75) is harder to compute since tt0 is not directly 
developed in harmonics. It can however be brought to a simpler form: 


(Ua, Vb[uU60]) = Vbu(Ma6,0) -|-VbU(1, 0) 

2 ^ ^ n 

= dd^bV^y-^Qab + VaV-yf 
d d +2 d 

212 ^ ^ ^12 


(77) 
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Note that V always applies to everything on its right, unless specified otherwise. One 
then obtain for the evoluation of p 

Pa = -'^a[vp\ - {Dr{d - 1) + a)pa + V • [AVpa] “ (78) 

Next, one projects ( [6^ on Mab to obtain 

{Mab,i^) = ^^Qab (79) 

a a + 2 

for the left hand side, and for the r.h.s. 


- {Mab,'^VU'tjj) + Dr{Mab,^u1p) + V • [AV(Mab,'i/')] “ a{Mab,i’) (80) 


Again, the last two terms are simple since {Mab,'ip) = 2f2Qab/[d{d + 2)]. Only the 
second angular/spherical harmonics survive the projection on Mab', they are eigen¬ 
vectors of Z\u, with eigenvalues —^(£ -|- 1) = 2d for d=3 and —= — 2d for d = 2. 
Hence the rotational diffusion term contributes a factor 


17 2 

Dr{Mab,2\u1p) = 


(81) 


Again, the first term of (80) is slightly more difficult to handle and has to be split as 


-(Mab, Vc'UUcV') = cV{Mab,UaUdPd) - ^cV{Mab,Ua0[P>]) 

Q 
d 


= -VcVPd{Mab, {Mad + ^)) - VcU^Xabc (82) 


d 

— c'^Pd^abcd ^c^^Xabc 

a a 

where we have introduced ^Xabc = {Mab,Uc&[ip]) which comes from the projection 
of uOlip] on the second harmonics and whose precise expression we won’t need. 
Putting everything together then yields Eq. (481 of the main text 


Qab — 


d-b2 


^c'iiPdddabcd 2dDj-Qab ^ ’ [ddt^Qab] ab ^c^Xabc (^^) 


with X = (d-b2)x/2. 


B From microscopic to mesoscopic Ito-Langevin equations 

ltd drift with functional dependences. The total gradient of D with respect to r^, 
notated as D'{r^, [p]) for lack of a better notation, is given by the chain rule 

A(r^, [p]) ^ ViZ7(r^, [p]) + J dr'^^^Vr,p(r') (84) 

where the first term on the r.h.s. comes from the explicit dependence of D on 
and the second from the dependence of D on p(r') which itself depends on r^. By 
convention, ViZ7(r^,[p]) is thus a “partial gradient” which acts upon the explicit 
dependence of D on (its first argument) but not on its implicit dependence through 
p. Using the explicit expression of p(r'), one has 

Vr,,p(r') = “ "'j) = Vr/(r' - r^) = -VrA(r' - r^) (85) 

j 
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Equation (85) shows why the notation Vr^, which indicates that the gradient acts on 
and not on r' or r,-, is essential—though cumbersome—to get the correct result. 
Inserting (85) into (84) and integrating by parts, one then finds that 


[p]) = Vi£)(r^, [p]) + J dr'(5(r'- r^)Vr'- 


[p]) 

Sp{r') 


= ViZ?(r^,[p])+ V 


SpiVf,) 


( 86 ) 


[p]) 


As a consistency check, and to understand how the operator 7p(F^ works, let 
us consider a simple example where the diffusivity of a particle at position is a 
linear functional of the density of its neighbours, convoluted by some kernel K: 

D{r^,,[p]) = J dr'p{r')K{r^-r') = J2K{r^-rj) (87) 

j 

In such a case, the total gradient of D with respect to can be directly computed: 


D'irf,, [p]) = - Tj), 


( 88 ) 


where VAT is the standard gradient of the function K[r). Applying the formula (86), 
one finds 


D'{r^,,[p\)= j dr'p{v')VK{Y^-v') + J dr'[Vr^^^^]K{r^ - r') 


(89) 


where the first term comes from the explicit derivative and the second one from the 
functional derivative. Note that, in the latter term, Vr^ does not act upon A'(r^ —r'). 
The computation can now be readily pursued, to give 

D'{r^,[p])= J dr'Y^6{r' -rj)WK{r^-r')+ J dr'[Vry (r'- r^)]A:(r^ - r') (90) 


Using again that Vr i5(r' — r^) = —Vr'5(r' — r^) and integrating by parts, one gets 


[p]) = ^ VK(r^ - r,) - Vif (0) = ^ V7f(r^ - r,) 
j jAm 


(91) 


Interestingly, the functional derivative generates a term — VAr(O) which will be 
absent whenever K{r) is a symmetric kernel or when the particles are not self¬ 
interacting. (By this we mean that the diffusivity D of particle p is a function of 
and a functional of p(r) — (5(r — r^); as such it is functionally dependent only 
on the density field of other particles, which differs from the total density by the <5- 
function self-term.) For complete generality however, we must retain this term in the 
Langevin equation (55) and throughout the construction of the Langevin equation for 
p(r), in order to derive a Fokker-Planck equation which is properly ordered as this is 
crucial to get the correct steady-state. Since the contribution of this term vanishes in 
many cases, it is often silently omitted in the literature. 

Let us further note that, looking at the simple derivation of equation (88), all 
ambiguities can be overcome by coming back to the microscopic definition of p(r). 
As often with functionals, the notational problem is present only when working at 
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the field level, where it stems from an underlying ambiguity concerning the gradient 
symbol V, which is used in different contexts to represent either a total derivative 
(acting on both implicit and explicit dependences of D{r, [p]), and hence containing 
in practice a functional derivative) or a partial derivative acting only on the first 
argument of D. This ambiguity translates into an ordering problem at the level of the 
Fokker-Planck equation. However, at the large deviation level—which corresponds 
here to the large size limit—we are effectively interested in the small noise limit 
of our stochastic partial differential equations. The various ordering of the Fokker- 
Planck equation differ by terms, similar to those that distinguish ltd and Stratonovich 
time-discretisation, which are generally negligible [51]. This is why the issue was for 
instance (rightly) neglected in |7I] . 


Mesoscopic Langevin equation. Let us now consider a function /(r) and compute 
the time evolution of /(r^(t)) where r^(t) is solution of the Langevin equation (55). 
Using Ito’s formula [Tgrfi] . one finds 


= [A + V^r]]\/f{rf,) + DAf{rf,) 


where 


A = V + ViZ?(r^,[p]) + (V, 




)-D(r^, [p]) 


(92) 


(93) 


Introducing = 5(r —r^) and using that, for any function H, H{y^) = f drp^iJ(r), 
Eq. (92) can be rewritten 


= J drp^(r,t)[(A-|-VTDr 7 )V/(r)-bi:)(r, [p])Z\/(r)] 


(94) 


where the derivatives in V/(r) and k\/(r) are now taken with respect to r. The 
last term in (94), even though it looks harmless, requires some explanations. Indeed, 
D(r^, [p]) hasaependencies on both through its explicit dependence and through 
its functional dependence on p(r') = ~ '^a)- When going from (92) to (94), 

we use that 


(5(r - r^)I?(r^, [p])Z\/(r^) = (5(r - r^)i:>(r, [p])Z\/(r), 


(95) 


where we have replaced all the r^’s by r but the one in p, i.e., p(r') = ~ ’’a) 

has not been replaced by 6{r' — r) -|- ~ !■«)• There is thus no dependence 

of L) on r through p. From now on, gradients VD{r, [p]) are not ambiguous anymore: 
they solely apply to the first argument of D which is the only place where D depe nds 
on r. There is thus no need anymore for the notation Vi. Integrating by part Eq. (94) 
then leads to 


= J dr/(r)V- [-{A + VWr])pf,{r,t) + W{Dpf,)] 

Alternatively, /(r^) can also be written 


/(r^) = / drp^(r)/(r) 


Since the equations (96) and (97) hold for any function /, one gets 


(96) 


(97) 


Pn = '^[-{A + V^'n)p^{r,t) +ViDp^)] 


(98) 












Will be inserted by the editor 


29 


Introducing the Gaussian white noise 

y/2DpA = E 

whose statistics satisfy 

(.1) = 0; (^(r, t).l(r', t')) = 5{t - t')5{v - r') 

one gets 

p(r) = V[-p(r)V + D^p{v) - P{^rj^)D - y^vl] 
which is Eq. (56) as required. 


(99) 


( 100 ) 

( 101 ) 


C Functional Fokker-Planck Equation 


The easiest way to derive the Fokker-Planck equation (58) starting from the Langevin 


equation (56) is to spatially discretize the latter in one spatial dimension 


Pi — Ai + -(ilj_i_|_i77i_|_i + 


( 102 ) 


where 


. ^ r T. 1 , 1, dDi dDi ■ 

A, - V. [ - + -D,{p^+, - p.-i) + -(^ - 


Pi\ 


B. 


1 


i,i+l — ~2'/‘^^i+^Pi+^ 


and Bii_i — ^ \^‘ 2 Dj^_\p 2 —\ 


(103) 

(104) 


Note that we use centred differences ViOi = |(Oi+i —Oi_i) to respect the isotropy of 
the equation. One can then use the general relation between ltd-Langevin dynamics 
and the Fokker-Planck equation m. 

Xi = A^ + B^jPj P{x) = — 2 E (105) 

i * i,j,k * ^ 


to derive the Fokker-Planck equation satisfied by P{pi). To do so, we first note that 


2 


BikBjk — Di+lPi+lSk,i+l + \/DiZiPiZi6k,i-l 

k k 


X \/ Dj+ipj+i5k,j+i + a/ Dj-ipj-iSk,j-i 
— ^ Pi+lPi+lPj + lPj + l^i,j ^Di — lPi—lPj + lPj + l^i—2^j 

+ ^/D, -iPi-i2Dj-ipj-i5ij 


(106) 

(107) 


To (slightly) lighten the notation, we write dp^ for the operator ^. The second order 
differential operator then becomes 


2 ^ ( dp^dp^Bij^B ^ ( Qp^ \^piPi+lPi+l ^pi — ’^Pi—lPi—l 

i,j,k i 

^ ^ ^pi {_^pi ^pi+ 2 ) ^ ^ ^pi (^pi-2 Di—ipi—i 

i i 
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Shifting by ±1 the indices in the two sums, one finally gets 


ij,k i 


( ^ 

9 ] 

( ^ 

9 ] 

\dp,-i 

dpi+i J 

\dpi-i 

dpi+i J 


DiPi 


(108) 


We now have to compute the drift term in the Fokker-Planck equ ation dp^Ai. 
Noting that Ai is a gradient, Ai = where Ci can be read in (103), the drift 

term can be rewritten 


2Y,dpA. = Y.9pM+i - C.-i) = 

i i i 

_ Pi \ ( Win ^*+1 ~ ^*-1 I ^ 

— / \ dpi--i ( piVi + Di + 


(109) 


2 \dp. 


'i+l 


dP^ 

dpi-i 


Putting everything together yields 


r) — F) 

P _ C'pi-1 Wi + l 


Pi / dP 

PiVi — PiViPi —— ' 


dP, 


2 \dp^+i dpi-i 


d - d 

^Pi — l ^Pi^l 


Pi Pi 


( 110 ) 


Let us now note that pi commutes with 9p._i — but not Pi, since the latter is 
in principle a function of all the Pj's. However, the application of dp^_-^ — 9pi^.i on Pi 
exactly cancels the one before last term in (110), so that 


^ = E - ^p.+i) 


p 


Taking at this stage the continuum limit and integrating once by part yields the 
correct Fokker-Planck equation 


P = — da;- 


d 


dp{x) 

where one has recognised the operator 


V p{x)V{x, [p]) - P{x, [p])Vp - P{x, [p])pVa 


Sp{x) 




'Sp{x) 


P (111) 


( 112 ) 


Equation (1121 is actually the best way to make sense of this operator and to see that 
the gradient really applies to the field with respect to which we are taking a funct ional 
derivative and not to the argument of this functional derivative. The equation (111) 
directly generalizes to (58) in higher dimensions. We note that to get the correct 
ordering between pP and Vr Fokker-Planck equation, it was necessary to 

correctly account for the atypical form of the spurious ltd drift as defined in the main 
text. 
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